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Abstract 

We determine the strong-isospin violating component of the neutron-proton mass difference from 
fully-dynamical lattice QCD and partially-quenched QCD calculations of the nucleon mass, con- 
strained by partially-quenched chiral perturbation theory at one-loop level. The lattice calcula- 
tions were performed with domain-wall valence quarks on MILC lattices with rooted staggered 
sea-quarks at a lattice spacing of b = 0.125 fm, lattice spatial size of L = 2.5 fm and pion masses 
ranging from m n ~ 290 MeV to ~ 350 MeV. At the physical value of the pion mass, we predict 
M n - M p \ a ~ U = 2.26 ± 0.57 ± 0.42 ± 0.10 MeV where the first error is statistical, the second error 
is due to the uncertainty in the ratio of light-quark masses, n = m u /md, determined by MILC [lj], 
and the third error is an estimate of the systematic due to chiral extrapolation. 
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I. INTRODUCTION 



It is a basic property of our universe that the neutron is slightly more massive than the 
proton. The electroweak interactions are responsible for this mass difference, which receives 
contributions from two sources. The strong isospin breaking contribution (also known as 
charge-symmetry breaking, for a review see Ref. j2j) is due to the difference in the masses of 
the up and down quarks, ultimately determined by the values of the Yukawa couplings in the 
Standard Model of electroweak interactions and the vacuum expectation value of the Higgs 
field. The other contribution arises from the fact that the proton and neutron carry different 
electromagnetic charges. The experimental neutron-proton mass difference of M n — M p = 
1.2933317 ± 0.0000005 MeV [j| receives an estimated electromagnetic contribution of [H 
M n — M p \ em = —0.76 ± 0.30 MeV, and the remaining mass difference is due to a strong 
isospin breaking contribution of M n — M p \ d ~ u = 2.05 =F 0.30 MeV. 

The recent suggestion of a vast landscape of possible universes emerging from string 
theory 0,0], and the previously hypothesized multi-universe scenario in general, have led to 
a resurgence of interest in the anthropic principle, and how it may provide a useful way of 
constraining fundamental parameters of nature. For this reason and others it is of interest to 
understand how quantities that affect nuclear physics and the production of elements depend 
upon the fundamental parameters: the length scale of the strong interactions, Aq C d, the 
light-quark masses, m u , m d and m s , and the electromagnetic coupling, a e . For instance, 
the basic fact that the neutron-proton mass difference is larger than the electron mass is 
central to the evolution and content of our universe (for a recent detailed discussion of the 
impact of an unstable proton see Ref. Q.). While a universe with a stable neutron would 
have a rich periodic table of nuclei, it is far from clear that the organic chemistry required 
for carbon-based life would exist @, H, |^ . 

In this work we have performed lattice calculations of the nucleon mass with domain-wall 
valence quark masses tuned to the staggered sea-quark masses of the MILC configurations 
(a mixed-action "QCD" calculation), and with valence quark masses that differ from the 
sea-quark masses (a mixed-action partially-quenched, "PQ", calculation). As the computer 
resources do not presently exist to perform such calculations at the physical values of the 
light-quark masses, they are performed at QCD pion masses of m^ tt ' ~ 290 MeV and 
m iatt. ^ M e y anc [ a ran g e f partially-quenched pion masses in between. These re- 
sults are combined with the formal framework of heavy-baryon chiral perturbation the- 



ory fHBvPT llO. Ill]) and partially-quenched heavy-baryon chiral perturbation theory 
(pqhbxpt EMM! E3), and the recent precise lattice extraction of the light-quark 
mass ratio, m u /md = 0.43 ± 0.01 ± 0.08 1] by the MILC collaboration, to calculate the 
strong isospin-breaking contribution to the neutron-proton mass-difference. 



II. THE FORMAL FRAMEWORK 



Nucleon observables have a systematic loop expansion about the limit of vanishing quark 
masses and external momenta defined by HB^PT [Tol 11 1. Extensive phenomenology has 
been performed with HB%PT, where it is found that some observables appear to converge 
quite rapidly in the chiral expansion, while others are less convergent. This state of affairs is 
magnified when dealing with three light flavors rather than two. Recently, the heavy-baryon 
formalism has been extended to describe the situation where the sea-quark masses differ 
from the valence quark masses [l2[ EH, EH HI; as is the case in partially-quenched lattice 
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calculations. 

The mass of the proton, which has valence quark content uud, when computed on con- 
figurations with sea-quarks, j and I, of mass rrij and mi has the form [15[ (we have used the 
opposite sign convention for the mass counterterms) 
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M p = M + (m u + m d ) (a + (3) + 2a (uij + mi) + - (2a - (3) (m u - m d ) 
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9an 



ml 



5m 3 ud + 3m% + 2m% + 3mf d + 2mf u + 3G VutVu + 6G VujVd + 3G Vdm 



m ]u + m L - m ld + 2m L + 3Gr, um + 3G VuiVu 



+ -g^r [ 5F ud + F uu + F ju + Fiu + 2F jd + 2F id + 2E mtVd + 2E Vu:Vu — AE Vu ^ Vd ] 1 , i 



where a,j3,a are the counterterms that enter at 0(m q ), qa is the nucleon axial coupling 
constant, g an is the A- nucleon coupling constant, and g± is the nucleon coupling to the 
flavor- singlet meson field, which makes no contribution in the QCD limit where mj —> m u 
and mi — * m d . The notation is such that m a b is the mass of the pseudo-Goldstone boson com- 
posed of quarks a and b. The functions G VatVb , and E VajVb are G VatVb = % Va rt b ( m i] a , m r) b , m x) 
and E, 



Va,Vb 



HvaVb( F vw F rib' F x)i respectively, where the function H Vanb is given by 
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and where the mass, m^, is given by m\ = 
Ft, = F(m ir , A, fi)) is given by 

F(m,A,/i) = (m 2 -A 2 ) ^VA 2 - m 2 log 



1 2a, m 

-m A log — 

2 \fi 2 



^(m^ + m 2 ). 



m^ + mfA. The A-loop function (e.g. 



'A - VA 2 -m 2 + ie 
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The nucleon mass has been computed to one higher order in the partially-quenched chiral 
expansion, 0(m 2 q ), inRef. 0, but with the small number of different quark masses available 
to us in our lattice calculation we will be unable to make use of that work. In the QCD 
limit, the proton mass computed from eq. (JTJ) becomes the QCD proton mass computed 
previously with HBxPT, 



M p = M + {a + [3 + 2a) {m u + m d ) + - (2a - (3) {m u - m d ) 



7>9a< + 



^9an 
3tt 



(4) 
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as required. The expansion of the neutron mass can be recovered from the expansion of the 
proton mass by interchanging the up and down quark masses, u <-> d. At the order to which 
we are working it is most convenient to replace the explicit quark masses in the expression 
for the proton mass with the leading-order expression for the pion mass to yield 



M +(a + (3 + 2a)ml - ~ (25 ~ P) (\ V 



8nf 



3 2 3 ^9 AN 77i 

- 2 9 A m, + 



where rj = m u jmd. The neutron mass is recovered by making the replacement r\ 
consequently 



M n - M p \ d ~ u = | (25-/3) Cj 



+ T] 



V \ 2 

mi 



(5) 

I/77, and 
(6) 



The one-loop contributions at 0(m^ 2 ) cancel in the mass-difference, as the pions are de- 
generate up to 0{m 2 q ). Analogous expressions for the partially-quenched proton masses can 
be found in Appendix 1X1 

The practical message one should take from the functional form of the proton mass in 
eq. (0) is that partial-quenching allows an extraction of isospin-violating quantities from 
isospin-symmetric lattices. 



III. DETAILS OF THE LATTICE CALCULATION AND ANALYSIS 

Our computation uses the mixed-action lattice QCD scheme developed by LHPC [13, 3 us- 
ing domain- wall valence quarks from a smeared-source on Nf = 2+1 asq t ad- improved 19f. l2p| 
MILC configurations generated with rooted 1 staggered sea quarks [26j] that are hypercubic- 
smeared (HYP-smeared) 0, H HI HJ . In the generation of the MILC configurations, the 
strange-quark mass was fixed near its physical value, bm s = 0.050, (where b = 0.125 fm 
is the lattice spacing) determined by the mass of hadrons containing strange quarks. The 
two light quarks in the configurations are degenerate (isospin-symmetric) . The domain- wall 
height is m = 1.7 and the extent of the extra dimension is L5 = 16. The MILC lattices were 
"chopped" using a Dirichlet boundary condition from 64 to 32 time-slices to save time in 
propagator generation. In order to extract the terms in the mass expansion, we computed 
a number of sets of propagators corresponding to different valence quark masses, as shown 
in Table I. On 468 bmi = 0.007 (denoted by Vi) lattices we have computed three sets cor- 
responding to the QCD point with a valence-quark mass of bm^wf = 0.0081 (Vi), three sets 
on 367 bmi = 0.007 lattices with a valence quark mass of bm dw f = 0.0138 (denoted by V 2 ), 
and two sets with a valence quark mass bm dw f = 0.0100 (denoted by V 3 ). On 658 of the 
bmi — 0.010 (V 2 ) lattices we have computed three sets at the QCD point with a valence- 
quark mass of bmdwf = 0.0138 (V2) and one set with a valence quark mass of bm^wf = 0.0081 
(Vi). The parameters used to generate the QCD-point light-quark propagators have been 



1 For recent discussions of the "legality" of the mixed-action and rooting procedures, see Ref. |2lL I22L \2c 
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Ensemble 


Theory 


bmi 


bm s 


bm dw f 


10 3 x bm res 


# props. 


2064f21b679m007m050 


QCD 


0.007 (VI) 


0.050 


0.0081 (Vi) 


1.604 ±0.038 


468x3 


2064f21b679m007m050 


PQQCD 


0.007 (Vi) 


0.050 


0.0138 (V 2 ) 


1.604 ±0.038 


367x3 


2064f21b679m007m050 


PQQCD 


0.007 (Vi) 


0.050 


0.0100 (V 3 ) 


1.604 ± 0.038 


367x2 


2064f21b679m010m050 


QCD 


0.010 (V 2 ) 


0.050 


0.0138 (V 2 ) 


1.552 ±0.027 


658x3 


2064f21b679m010m050 


PQQCD 


0.010 (V 2 ) 


0.050 


0.0081 (Vi) 


1.552 ±0.027 


658x1 



TABLE I: The parameters of the MILC gauge configurations and domain-wall propagators used 
in this work. For each propagator the extent of the fifth dimension is L5 = 16. The notation of 
quarks, V%, V2, V3, is defined in the text. The last column is the number of propagators generated, 
and corresponds to the number of lattices times the number of different locations of sources on 
each lattice. 

"matched" to those used to generate the MILC configurations so that the mass of the pion 
computed with the domain- wall propagators is equal (to few-percent precision) to that of the 
lightest staggered pion computed with the same parameters as the gauge configurations jib^ . 
The lattice calculations were performed with the Chroma software suite (3ll . l32j ] on the high- 
performance computing systems at the Jefferson Laboratory (JLab). Various proton and 
pion correlation functions were constructed from the three (distinct-mass) light-quark prop- 
agators that were generated on the bm\ = 0.007 lattices, and the two light-quark propagators 
that were generated on the bm\ = 0.010 lattices. Differences between the various proton 
states were also constructed. It is useful to define the proton mass splitting: 

AM p (V a , V b , V c ; V d ) = M p (V a , V b , V c ; V d ) - M p (V d , V d , V d ; V d ) , (7) 

where the indices a, b, c range over 1, 2, 3. Effective mass plots for the ratios of proton 
correlators that give rise to these mass splittings are displayed in figs. [TJ |2] and 01 Results 
for the extracted proton mass splittings are given in Table |H] and displayed in fig. 0] The 
various pion masses relevant to the analysis are also shown in Table El Using this data in 
conjunction with the PQHBxPT formulas given above and in the appendix, the coefficients 
<5> @, gAN, and g\ that appear in eq. (0) were extracted using various nonlinear fitting 
techniques. A tree- level, 0(m q ), analysis of the mass differences shown in Table ITD using 
the expressions for the partially-quenched proton masses given in the appendix, allows for 
an extraction of the isospin-breaking coefficient (2a — (3j /3, and the isospin-conserving 

coefficient a + (3, as shown in Table IIHI This leads to a prediction for the strong isospin 
breaking at the physical quark masses, at 0(m q ), of 

M n -M p \ d ~ u = 1.96 ± 0.92 ± 0.37 MeV , (8) 

as shown in Table IIV1 The first error is statistical and the second error is due to the 
uncertainty in the determination of the ratio of light quark masses, rj = m u /m d = 0.43 ± 
0.01 ±0.08 by MILC Q. 

The analysis of the partially-quenched data at 0(m^ 2 ) introduces contributions to the 
various proton masses that vanish in the QCD limit, a characteristic feature of the partially- 
quenched theory. In this case there is a contribution from loop diagrams involving the 
flavor- singlet field — due to a mismatch between the ordinary flavor singlet and that of 
the graded-group — which depends on an axial coupling constant, g\ [l5|, that must also 



5 



0.06 - 
0.04 - 
0.02 - 



♦ * * i i 



10 



12 



0.06 - 
0.04 - 
0.02 - 



• ■ * m 5 :: 



10 



12 




t (l.U.) 

FIG. 1: Effective mass plots for the ratios of proton correlation functions that give rise to the mass 
splittings AM p (V a , V c ; V\) with a, b, c=l, 2. 



be determined from the lattice data. Fortunately, with multiple partially-quenched proton 
states, gi can be extracted simultaneously with (2a — pj /3, a + /3, and gAN- However, 
given that a and (3 enter at tree-level, while g\ and gAN enter at the one-loop level, the 
fractional uncertainty in g\ and gAN will be parametrically larger than that of (2a — (3^ /3 

and a + j3 (assuming natural sizes). 

At this order there are contributions from loops diagrams involving nucleons and diagrams 
involving A's. The axial coupling between the nucleons and the pions is taken from the 
chiral perturbation fit to the recent LHPC lattice calculation [33|. We use a value of gA 
that is the average of = 1.25, corresponding to the QCD point on the bmi = 0.007 
lattices, and g^ = 1.24, corresponding to the QCD point on the bmi = 0.010 lattices. 
Similarly, for the pion decay constant, we use a value of f n = 149.8, which is the average 
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FIG. 2: Effective mass plots for the ratios of proton correlation functions that give rise to the mass 
splittings AM p (V a , V c ; Vi)with a, b, c=l, 3. 



of f n = 147.8 MeV corresponding to the QCD point on the bmi = 0.007 lattices, and 
f n = 151.8 MeV corresponding to the QCD point on the bmi = 0.010 lattices. We fit <?atv to 
the data, and find \g^ N \ = 0.60 ± 0.66. This value is slightly smaller than the SU(6) value 
of gAN — ~%9a 2 ■ The mass difference between the A and nucleon must be input into the 
extraction. Ideally, this would also be fit to the data, but we do not have precise enough 
data to permit such a fit. As the neutron-proton mass difference is highly insensitive to this 



2 A three-parameter fit with |<7ajv| = 1.8 (the value obtained at tree-level from A — > Ntt) yields M n — 
M p \ d ~ u = 2.06 ± 0.99 ± 0.38 MeV while fitting with nucleon loops alone (g A N = 0) gives M n - M p \ d ~ u = 
2.29 ±0.53 ±0.43 MeV. 
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FIG. 3: Effective mass plots for the ratios of proton correlation functions that give rise to the mass 
splittings AM p (V a , V c ; V2) with a, b, c=l, 2. 

quantity, we use the experimental value Ma — Mn = 293 MeV. At one-loop order, 0(m^ 2 ), 
we find 

M n — M p \ d ~ u = 2.26 ± 0.57 ± 0.42 ± 0.10 MeV , (9) 

where the last error is an estimate of the systematic error due to truncation of the chiral 
expansion. It is reassuring that the predicted neutron-proton mass difference is relatively 
insensitive to the order in the chiral expansion, as shown in Table IIVI Both the tree-level 
and the one-loop extraction of the neutron-proton mass differences are consistent with the 
"experimental" value of M n - M p \ d ~ u = 2.05 ± 0.30 MeV. 

An interesting observation can be made by comparing the proton mass differences on 
the two different lattice sets, as shown in Table HI and displayed in fig. HI Within errors, 



a -0.02 
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-0.04 

n 



s 



Quantity 


Mass (Difference) (I.u.j 


Mass (Difference) (MeV) 


Fitting Range 


m^Vi, Vi; Vi) 


0.1864 ± 0.0011 


294.2 ± 1.7 


5 — ► 15 


m^iVx.V^Vx) 


0.2066 ± 0.0010 


326.2 ± 1.6 


5 — > 15 


m n (V2,V 2 ;Vi) 


0.22473 ± 0.00091 


354.4 ± 1.4 


5 —> 15 


^3; vi) 


U.lyzy ± U.UUlz 


304.5 ± 1.9 


5 — > 15 


m*(V 3 ,V 3 ;Vi) 
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n 1 s/m -i- n nm ^ 
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zyi.u m z. i 
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— ► 10 


m n (V 1: V 2 ;V2) 


n 9n^n -i- n nm 9 

U.ZUOU HZ U.UU1Z 


QOQ 7 _|_ 1 n 
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r. , it; 
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n 99^^ -i- n nm 1 

U.ZZOU ZIZ U.UUII 


QfO Q _|_ 1 O 

ooz.y in 1.0 


K i it; 

— > 10 


AM p (Vi, V\, V 2 


Vi) 


0.0163 ± 0.0019 


25.7 ± 3.0 


5 — » 12 


AM p (V 2 ,V 2 ,Vi 


VI) 


0.0209 ± 0.0029 


32.9 ± 4.7 


5 — » 12 


AM P (V 2 ,V 2 ,V 2 


VI) 


0.0353 ± 0.0041 


55.8 ± 6.5 


5 — > 12 


AM p (Vi,Vi,V3 


VI) 


0.0049 ± 0.0010 


7.7 ± 1.6 


5^11 


AAf p (v 3 ,v 3 ,vi 


Vi) 


0.0061 ± 0.0016 


9.7 ± 2.5 


5 — > 11 


AAf p (v 3 ,v 3 ,v 3 


Vi) 


0.0109 ± 0.0024 


17.2 ±3.8 


5 -» 11 


AM p (Vi, V\, V\ 


v 2 ) 


-0.0309 ± 0.0038 


-48.8 ± 6.0 


4^ 11 


AM p (Vx, Vi,V 2 


v 2 ) 


-0.0161 ± 0.0022 


-25.5 ±3.5 


4^ 11 


AA/ P (v 2 ,v 2 ,vi 


v 2 ) 


-0.0137 ±0.0016 


-21.6 ±2.6 


5 12 



TABLE II: The pion masses and proton mass differences calculated on the bmi = 0.007 and 
bm\ = 0.010 MILC lattices. The notation of valence and sea quarks, Vi 5 2,3, is defined in the text. 
A lattice spacing of b = 0.125 fm has been used. 



Extraction 


l(2a-(3) (l.u.) 


a + ]3 (l.u.) 


9i 


\9an\ 


X 2 /dof 


LO 0{m q ) 


0.198 ± 0.093 


2.07 ±0.08 






0.56 


NLO 0{mf 2 ) 


0.229 ± 0.058 


3.4 ± 1.1 


-0.10 ±0.35 


0.60 ± 0.66 


0.21 



TABLE III: Parameter Table. The values of the parameters in the partially-quenched chiral La- 
grangian as determined by a ^-minimization fit of the theoretical proton mass differences given 
in Appendix lAl to the lattice data given in Table ILT1 The isospin-conserving combination of coun- 
terterms, a + f3, is renormalization-scale dependent. We have renormalized at fj, = 1 GeV. 



the magnitude of the mass differences are independent of the value of the sea-quark mass. 
This is consistent with the leading-order chiral expansions given eq. (0) and in Appendix 1X1 
Higher-order contributions to these mass differences in the chiral expansion, which give rise 
to deviations from these equalities, will be become more visible with increased statistics. 

There will be finite lattice spacing contributions to the parameters that we have extracted 
in this work. The recent developments in the inclusion of finite-lattice spacing effects in 
mixed-action theories in allow us to determine where such corrections enter and to 
estimate how big the corrections should be. The lattice spacing is introduced into the 
mixed-action theory by extending the SU{2)l ® SU{2)r lie-algebra to a graded lie-algebra 
that makes the distinction between sea and valence quarks explicit. The lattice spacing is 
incorporated by a spurion field with the appropriate transformation properties under the 
graded group, e.g. see Ret. HHiH l37| . There is a leading-order contribution at 0(a 2 mj) 



to the nucleon mass (where we are assuming that the exponentially-suppressed contribution 
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FIG. 4: The 'partially- quenched proton mass differences (in MeV) calculated from the bm\ = 0.007 
and 0.010 MILC lattices plotted vs the pion mass composed of sea quarks. Various data have been 
displaced horizontally by small amounts for display purposes. A lattice spacing ofb = 0.125 fm has 
been used. 



Extraction 


M n - M p \ d ~ u (MeV) at m^ 8 - 


LO 0{m q ) 


1.96 ±0.92 ±0.37 


NLO 0(m q 12 ) 


2.26 ±0.57 ±0.42 



TABLE IV: The neutron-proton mass-splitting at the physical value of the pion mass, mP hys ' = 
140 MeV, extracted from this partially-quenched lattice calculation, using the parameters shown 
in Table II I II The lattice spacing used to convert between lattice units and physical units is b = 
0.125 fm. The first error is statistical while the second error is due to the uncertainty in the ratio 
of quark masses, m u /md, in the MILC calculation 0. 



at 0(a mjj) from the finite L5 is numerically insignificant). However, such terms do not 
contribute to the mass differences between the proton states that we have used to extract 
the parameters. Finite lattice spacing contributions to the nucleon mass that depend upon 
the light-quark masses start at 0(a 2 m q ). Therefore, we expect the finite lattice spacing 
corrections to our results to be parametrically suppressed and small. In contrast, the finite 
lattice spacing contribution to the nucleon mass itself is expected to be roughly the same 
size as the contribution from the <r-term, rendering an extraction of the nucleon a-term 
(e.g. see Ref. [38]) somewhat unclean if the physical value of the nucleon mass is used 
in the extraction. Finite-volume effects in the baryon mass splittings are estimated to be 
negligible [39^ . 
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IV. CONCLUSIONS 



In this work we have performed the first lattice calculation of the neutron-proton mass- 
difference arising from the difference between the mass of the up and down quarks, and 
find M n - M p \ d ~ u = 2.26 ± 0.57 ± 0.42 ± 0.10 MeV. This value is consistent with the 
number 2.05 ± 0.30 MeV based upon the experimentally-measured mass-difference and the 
best estimate of the electromagnetic contribution. It is clear that further lattice calculations 
are warranted in order to make a precise prediction for this quantity, which will then enable 
a precise determination of the electromagnetic contribution to this mass- difference 3 . 

The neutron-proton mass difference is but one of the manifestations of charge-symmetry 
breaking that has been a focus of both theoretical and experimental investigations for many 
years |2j. With the recent progress in extracting the nucleon-nucleon scattering amplitudes 



from fully-dynamical lattice QCD |41[, one can imagine using partially-quenched calculations 
to extract the charge-symmetry breaking contribution to nucleon-nucleon scattering, and 
other processes, in future investigations. 
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APPENDIX A: NUCLEON MASS IN PQHBxPT 

In this appendix we give the expression for a proton composed of valence quarks 
V\ViVi, ViViV 2 , V 2 V 2 Vi and V 2 V 2 V 2 on an isosp in- symmetric sea composed of two light 
quarks V 1: V\. The proton masses are, starting with the QCD point, 

~ m Vi,Vi;Vi 



M p (Vi, Vi, Vi, Vi) = M + (a + P + 2a 
_^9a_ 3 9an tp 

~ i5^2 m Vi,Vi;Vi - g^2p^Vl,Vl;Vl 

M p (Vi, V x , V 2 ; Vy) = M + 2 a m 2 Vl y. Vl + ^ (5a + 20) m~ VuVl . Vl 



2 



+ l(a + 4/3) 



g \ - / m V- 2 ,V2;Vi 



ml 



24nf 2 \2 Vl ' Vl ' Vl ""VuVaVi 
9a9\ {5 3 3 3 3 



384tt/ ; 



9an 



72tT 2 / 



(lAm\ yi . Vi + 4:m 3 VuV2 . Vl - 27m^ 2)V2 . Vl + 9m^ yi . F m^y 2 . Vl ) 
( 2i 7 Vi,yi;Vi + 9Fv u v 2 ;Vi + Fv 2 ,v 2 ;Vi + m v 1 y 1 ;Vi S y2,v 2 -y 1 - m v 2 y 2 y 1 Sv 2 y 2 y 1 ) ; 



M P (V 2 , V 2 , Vi, Vi) = M + 2am 2 VuVl . Vl + ^ (557 + 2(3) 



g V " ; m V 2 ,V2;Vi 



+ 6 ( s + 4/9 ) m ^;^ 



96tt/ 2 

96tt/ 2 
,2 



(20m^ iy2 . yi + H iyi . Fl m F2 ,y 2;Fl - llrr4 2>V2;Vl ) 
( 4m Fi,v 2 ;yi + ^.v^Vi^.VajV! - l ?>™v 2 y 2 y) 



9i ( 1 0^,3 



384tt/ 2 

9an 
72tt 2 / 2 



(l8m^ ljVL;Vi - 4:m^ uV2 . Vl - 23rr^ a y 2 . Vl + 9m 2 VuVi . v m V2 y 2 . Vl ) 

( S-^i.VisVi + 7F Vli y 2;yi + 2Fy 2) y 2; vi + ^y^ySVa.VaiVi - ^Va.VajVi'S'va^lVi) ! 



M p (y 2 , F 2 , V 2 ; VI) = M + 2am Vl)yi . yi + (« + /?) m 



V2,V 2 ;Vi 



,2 



96tt/ 2 



iSm^y^y + 9m 2 Viyi . Vi m V2 y 2 . Vl - lm\ 2 y 2 y 



9a_ 
P 

( Am yy 2 y + 9m yyy m v 2 yy - 13my 2) y 2;Vl ) 



48tt/ 



2 



2 



+ 9m yy 1 y 1 m v 2 y 2 y 1 ~ ^ m \ 2 y 2 y x 



96tt/ 2 

Fv 1 y 2 yi + F V2 y 2 . Vl ) , (Al) 



9an 
12tt 2 / 
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where we have used the leading-order relation between the quark masses, and the function 
S n = S(m n , A, /i) is 

'A - VA 2 -m 2 + ie\ . ( ( m 2 \ I 



S (m , A, „) = yA^.o g { A + ^_y +i( ) ~ A (>o g + 3 j . (A2) 

The function F is defined in the text. 

The meson masses at leading order in the chiral expansion, starting with the QCD point, 

are 

ml(Vi,Vi;Vi) = rn 2 Vi Vi . Vi = 2 A m Vl 
ml(V u V 2 ;V 1 ) = m 2 ViV2 . Vl = A (m Vl + m V2 ) 

ml(V 2 ,V 2 ;V 1 ) = m 2 V2>V2 . Vl = 2 A m V2 , (A3) 

where A is a strong interaction mass-scale, and the semicolon separates valence and sea 
quarks. 
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